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Abstract We present a detailed model of stray-light suppression in the spectrometer 
channels of the Ultraviolet Coronagmph Spectrometer (UVCS) on the SOHO space- 
craft. The control of diffracted and scattered stray light from the bright solar disk is one 
of the most important tasks of a coronagraph. We compute the fractions of light that 
diffract past the UVCS external occulter and non-specularly pass into the spectrometer 
slit. The diffracted component of the stray light depends on the finite aperture of the 
primary mirror and on its figure. The amount of non-specular scattering depends mainly 
on the micro-roughness of the mirror For reasonable choices of these quantities, the 
modeled stray-light fraction agrees well with measurements of stray light made both 
in the laboratory and during the UVCS mission. The models were constructed for the 
bright H I Lya emission line, but they are applicable to other spectral lines as well. 

Keywords: Instrumental Effects; Spectrum, Ultraviolet 
1. Introduction 

Until the 20th century, total solar eclipses were the only means of observing the hot solar 
corona. However, with the invention of the internally-occulted coronagraph by Bernard 
Lyot in the 1930s (Billings, 1966';'Koutch my, 1988] ), the addition of an external occulter 



( |Evans, 1948 ; Newkirk and Bohlin, 1963 1, and the development of an ultraviolet coro- 
nagraph spectrometer in the 1970s ( |Kohl et ah, 1978[ |Kohl et al, 1980| l, a continuous 
and detailed exploration of coronal plasma physics became possible. Such instruments, 
combined with spectroscopic diagnostic techniques, have become powerful tools for 
measuring a wide range of plasma properties in the acceleration regions of solar wind 
streams and coronal mass ejections (see, e.g., INewkirkI 119671 IWithbroe et al. \ 119821 
IKohlef aZ.i r2006; Ho ward ef a/ . 2008; C ranmerll2009] l. 

One of the most demanding requirements of a successful coronagraph is the suppres- 
sion of "stray light" scattered from the bright solar disk. The coronal emission tends 
to be many orders of magnitude less bright than the emission from the Sun's lower 
atmospheric layers. At a heliocentric distance of two solar radii (i?©), the extended 
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corona is approximately 10^^ times less bright than the disk at visible wavelengths. 
In the ultraviolet the corona is relatively brighter, but still is only about 10^^ (at H I 
Lya 1216 A) to 5 x 10^^ (at Mg X 610 A) times the disk intensity. These numbers 
correspond to the lowest-intensity regions off the solar limb (coronal holes). Because 
Ught rays coming from the disk are separated in angle by only fractions of a degree 
from the rays coming from the corona, there are many ways that a small fraction of the 
former can easily contaminate the latter. A key purpose of an efficient coronagraph is to 
suppress as much of this stray light as possible through the use of re-imaging, occulters, 
baffles, light traps, and also by making the mirrors as smooth and accurate as possible. 

This paper presents a detailed analysis of the stray light properties of the Ultraviolet 
Coronagraph Spectrometer (UVCS) instrument (Ko hl et al, \995\ Kohl et al, \991) 



onboard the Solar and He liospheric Observatory (SOHO) spacecraft ( Domingo et al, 1995H 



Fleck and Svestka, 1997| l. At present, there are only a few published measurements of 
the stray light properties of UVCS, and many observations have been taken without 
dedicated stray-light measurements. Thus, there is a need for a model that predicts 
the stray-light intensities at ultraviolet wavelengths. The dominant source of UVCS 
stray light has been shown to be diffracted light from the external occulter that is 
non-specularly scattered by the telescope mirror into the spectrometer entrance slits. 
IKohl et a7n(ll995 i) described other possible contributions to the stray light and explained 
how they are controlled. This paper primarily describes the non-specular scattering 
from the telescope mirror, which accounts for essentially all of the UVCS stray light 
background. 

The model developed for this paper applies mainly to the Lya spectrometer channel 
of UVCS, which observes the strong H I Lya 1216 A emission line and other spectral 
lines within a wavelength range of 11 00 - 1361 A. The model also is valid for the O VI 
channel, which is optimized to observe the O VI 1032, 1037 A doublet, as well as other 
spectral lines within a wavelength range of 937 - 1 126 A in first order. The redundant 
Lya path of this channel can be used to observe from 1 166 to 1272 A in first order as 
well. The stray-light properties of the UVCS white-light channel (WLC) were described 
bv lRomoli et a/.l(J99 3 ) and are not discussed here. 

The outline of this paper is as follows. Section 2 presents an overview of the relevant 
UVCS optical paths and occulting surfaces. Section 3 describes how much light from 
the solar disk is diffracted by the UVCS external occulter, and Section 4 describes 
how much of this light enters the spectrometer slit via finite-aperture diffraction and 
non-specular scattering. Section 5 summarizes how the contributions from a range of 
locations on the UVCS primary mirror are summed to determine the total stray light 
intensity. Section 6 presents the results of a calculation for H I Lya stray light as a func- 
tion of the observation height and compares the model predictions with measurements. 
Section 7 contains a brief summary of the paper and a discussion of the capabilities of 
future instrumentation. Finally, the Appendix describes how the theoretical rates of non- 
specular scattering were derived separately for the effects of mirror micro-roughness 
and finite-aperture diffraction. 



2. UVCS Optical Path Geometry 

The UVCS instrument contains three reflecting telescopes that feed two ultraviolet 
toric -grating spectrometers and one white-light polarimeter ( |Kohl et al, 1995| l. Unlike 
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Figure 1. (a) Scliematic illustration of the UVCS/SOHO optical path for coronal rays (i.e., light coming from 
an infinite distance), (b) Optical paths for rays diffracted from the tip of the external occulter. The horizontal 
distance along the mirror [x] increases from left to right. In both panels, the mirror is the light-gray surface 
at the bottom, the external occulter is the dark-gray surface at the top (at a distance /^cxo above the mirror), 
and the spectrometer slit is the light-gray split suii'ace at a distance / above the miiTor. The internal occulter 
is the dark-gray suii'ace sitting on top of the miiTor and covering up all x > Xmax- 

many other coronagraphs, the UVCS occuhing surfaces are linear, not ckcular, so that 
the incoming rays are well matched to the linear geometry of the spectrometer entrance 
slit. Rays from the bright solar disk are blocked by an external occulter and enter a 
sunlight trap, whereas rays from the extended corona are reflected by the primary mirror 
into the spectrometer slits. These slits are oriented in the direction tangent to the solar 
limb. The coronal images can be positioned on the slit centers in heliocentric radius [r] 
anywhere between about 1.4 and 10 Rq and rotated around the Sun in position angle. 
The slit length projected on the sky is 40 arcminutes, or approximately 2.5 Rq in the 
corona, and the slit width can be adjusted to optimize the desired spectral resolution 
and count rate. 

Figure [T^ illustrates the geometry for light coming from the extended corona, being 
reflected by the UVCS primary mirror, and entering the spectrometer slit. The geometry 
is the same for the Lya and O VI channels. This diagram is not intended to be an 
accurate schematic of the instrument, however. We illustrate the mirror as being flat, 
but we take account of the fact that it focuses coronal light on the slit. The spatial scales 
are exaggerated in order to show how the relevant angles and distances are defined 
below. For example, all angles of reflection and scattering are small (i.e., near normal 
incidence on the mirror), with sm9 « tan 9 w 6 applying in most cases. Because the 
analysis of this paper is mainly concerned with relative departures from the specularly 
reflected ray paths, it is unimportant which of the rays in Figure[T]is portrayed as exactly 
normal to the mirror 

The external occulter is shown as the upper-most horizontal surface in Figure[T] It is 
positioned so that all specular rays from heliocentric distances less than a given distance 
Tcxo miss the mirror entirely and are absorbed in the UVCS sunlight trap (see dashed 
blue rays). Thus, we can only observe the extended corona at radii robs > ''cxo (see 
solid red rays). For UVCS, Toxo = 1-2 Rq. In Figure [T] the horizontal coordinate x is 
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measured along the length of the muTor. The origin x = is defined as the left edge 
where the ray coming from rcxo that just barely misses the external occulter also just 
barely misses the mirror 

The observation angle 6'obs, which corresponds to the angular distance between the 
desired observation height and the rays that just miss the mirror, is given by 



7obs — Oq 



^obs 



(1) 



where 9q is the angle subtended by one solar radius on the sky. For the mean heliocen- 
tric distance of the SOHO spacecraft, 0q = 4.698 x 10^^ radians = 0.269° = 968.9". 
For UVCS, the external occulter is positioned at a height of Hcxo = 170 cm above 
the mirror The maximum possible region of the mirror that can be illuminated by the 
corona extends from x ~ to x — iJcxo tan^obs- We can define these two limiting 
positions as xi and X2, respectively. The area of the mirror exposed to coronal light can 
be reduced by moving the UVCS internal occulter (which sits on top of the mirror) to 
the left of its default position by a distance x-mt- The default ("zero over-occulting") 
position corresponds exactly to X2, and most UVCS observations have been made with 
the internal occulter positioned at Xi^t = 1.5 mm to the left of X2- Thus, the so-called 
"vignetting function" of UVCS is defined by the radially increasing mirror area that is 
filled by coronal rays. For a given observation height, the unvignetted mirror area spans 
the distances between Xnun = xi ~ and a;inax — X2 — Xint- 

It should also be noted that neither X2 nor Xmax can exceed the actual physical size of 
the mirror in the x-direction, which for UVCS is 7.2 cm. The corresponding maximum 
observation height at which the mirror is "filled" is robs — 10.3 Rq. 

Parallel rays from an infinite distance are reflected by the mirror and are focused at 
the slit, which is positioned at / = 75 cm above the mirror For the orientation shown 
in Figure[T^, the left edge of the slit is fixed in place (exactly at the focus point) and the 
right edge of the slit is the part that opens and closes to achieve a given slit width w. In 
this analysis, all coronal rays - from the coronal region of interest - that are reflected 
from mirror positions between Xmin and Xmax are assumed to pass through the slit. 

FigureflJ? shows the behavior of light that is diffracted from the edge of the external 
occulter. From the point of view of geometrical optics, these rays all appear to originate 
approximately at the tip of the external occulter itself {i.e., at a finite distance away from 
the mirror). Thus, these rays are not focused to a point at /, but instead are focused at a 
distance i?diff above the mirror, which is given by 

' ' ' (2) 



and iJdiff = 134.21 cm for UVCS. Note that the specular rays in Figure[T}5 are nearly 
all blocked from entering the slit. The only exception would be the right-most rays that 
strike the mirror at X2 = i/cxo tan ^obs- These are eliminated by just a small amount 
of over-occulting by the internal occulter. However, the diffracted beams also have 
substantial energy at non-specular angles, and these are discussed further in Section 
4. 
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Figure 2. Fractional reduction in iiradiance due to diffraction at the external occulter, shown as a function 
of dimensionless distance [?;] behind the occulter. Laboratory measurements for a straight edge (points) and 
a serrated edge (dashed curve) are compared with computed diffraction curves (solid curves). 



3. Diffraction from the External Occulter 



The external occulter is the first line of defense against stray-light contamination. In this 
section, we compute how much light is diffracted by the external occulter and falls onto 
the mirror We ignore the diffracted light that does not hit the mirror, but instead strikes 
other structural components of the instrument where it can undergo multiple reflections 
before possibly finding its way to the entrance slit. Specifically, we ignore the stray 
light produced by scattering off the mirror edge and the surface of the internal occulter 
These effects can be important at visible wavelengths, but they have been shown to 
be negligible in the ultraviolet for externally-occulted designs like UVCS (see, e.g., 
IRomoli et a/.lfT993llKohl et a/.lfT995l l. 

The problem of parallel rays diffracting around an ideal one-dimensional straight 
edge (i.e., a semi-infinite screen) is treated in standard textbooks ("e.g.. IBorn and Woifl 
119991 iHechtl I2002l l . In practice, it has been found that serrated, or "toothed" edges per- 
form better than ideal straight edges ( jNewkirk and Bohlin, 1963[ Koutchmy and Belmahdi, 1987H 
IVerroi et al, 2008] l. In the late 1980s, laboratory measurements were made at the Smith- 
sonian Astrophysical Observatory for both straight-edged and serrated occulters, in 
combination with the Spartan 201 UV coronagraph spectrometer (see IRomoli et al.\ 
I1993I I. The laboratory measurements were done with a light source that subtends the 
same solid angle as the solar disk. The serrated occulter used on UVCS has symmetric 
sawtooth-shaped notches {i.e., resembling a repeating triangle waveform) with a sepa- 
ration of the peaks, parallel to the occulter edge, of 229 /im. The depth of the peaks, 
perpendicular to the occulter edge, is 267 /im. 
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Figure 12] illustrates laboratory measurements and theoretical predictions for the re- 
duction in irradiance {i.e., radiative flux) seen as an observer moves progressively into 
the occulter's shadow. This quantity is shown specifically for the H I Lya wavelength, 
and the distance behind the occulter was converted into dimensionless Fresnel diffrac- 
tion coordinates 



where A is the wavelength of interest and xq is the mirror position that corresponds to 
rays from the near edge of the solar disk that barely miss the external occulter. Because 
these rays do not strike the mirror, xq is negative in our coordinate system, and it 
is given by — i/oxo tan(6'Q[rcxo/-R0 — !])■ Thus, all values of v that are relevant for 
UVCS are positive. We denote the relative irradiance quantity plotted in Figure |2] as 
Div). 

The uppermost curve in Figure |2] corresponds to parallel rays (i.e., from a point- 
source on the sky at infinite distance). This is given by 



Div) 



I -Civ) 



l-Siv) 



(4) 



where Civ) and S'(u) are the standard Fresnel integrals ( [Abramowitz and Stegun, 1972| l.[ 
When observing the Sun, this function must be convolved with the finite size of the 
solar disk on the sky, since the total diffracted light is a combination of beams be- 
tween the "near edge" of the solar disk (i.e., only 0.2 Rq behind the 7-oxo ray that just 
grazes the left edge of the mirror) and the "far edge" (up to 2.2 Rq behind the grazing 
ray). For the H l Lya emission line, we consider the solar disk to be uniformly bright 
( |Curdt et al, 2008) 1. The result of this numerical convolution (lower solid curve) agrees 
very well with the measured data for a straight-edge occulter (points) taken with a light 
source with the same finite solid angle as the Sun (see also' Romoh et al.\\1993i . 

It can also be seen from Figure |2] that the serrated edge produced a significantly 
lower amount of diffraction than the straight edge. The dashed curve corresponding 
to the serrated edge is denoted Dsiv). The laboratory measurements extended up to 
an occulting distance corresponding to v ^ 38. UVCS observations made at robs < 
2.74 Rq fall into the range of heights where the laboratory measurements can be applied 
via interpolation. However, UVCS observations at larger heights require some degree 
of extrapolation beyond the measured range. For example, the largest height of j'obs ~ 
10 Rq fills the mirror with a range of v values between about 5 and 220. After some 
experimentation, we found that a power-law extrapolation curve that scales as best 
reproduced and continued the slope of the measured data. I Romoh et al.\ ( 1 19931 1 used 
an exponential function which may become lower than the i;""^ curve at large values 
of V. Thus, our power-law extrapolation may give rise to a slight (but conservative) 
overestimate of the diffraction at large observation heights. 



4. Non-specular Rays Entering the Slit 

Each point on the mirror between a;,nin and Xmax acts as a "source" of solar-disk light 
that has been diffracted by the external occulter. These specular rays are focused at a 
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Figure 3. Non-specular geometry for diffracted rays originating at an arbitrary point [x] on tlie mirror surface 
and passing tlirougli the spectrometer slit. 



point -ffdiff above the mirror, and are blocked by the left-side slit edge. However, the 
light is not concentrated all along the ideal specular ray, but instead has a characteristic 
non-specular spread that subtends a broader range of solid angle. Thus, there are al- 
ways non-specular rays from the mirror that are able to pass through the slit. Figure [3] 
illustrates those rays, which are denoted by non-specular angles ranging between 6'min 
and 0max- For a given point [x], the values of these angles are given by 

0n.in = tan-1 (J^) ^ tan-1 ( (5) 

\X2 - X J \X2 - X J 



= tan-1 (J^^) - tan-i ) (6) 

\X2 - X J \X2+W-Xj 

where w is the slit width. Note that for the special case of x — X2, some specular rays 
pass through the slit since 6^in — and 6'i„ax = tan^^(w//). Of course, for any 
finite amount of internal over-occulting, these rays never get reflected at all, since we 
endeavor to keep all rays that reach the miiTor at a; < Xmax < 2^2 ■ 

The fraction of reflected light that passes non- specularly into the spectrometer slit is 
computed by summing the effects of two distinct physical processes: /) diffraction due 
to the finite size of the mirror, and //) scattering due to mirror micro-roughness. The 
Appendix gives a derivation of these two components. We note that this separation into 
two parts is an approximation of the full angular distribution of non-specular intensity, 
which for UVCS is not known with sufficient detail to specify as a single function. It is 
useful to specify these two effects separately so that different assumptions about their 
origin and magnitude can be made without each "contaminating" the other Taking their 
direct sum may result in an amount of non-specular scattering that is slightly larger than 
the actual combined effect of these two processes, but this errs on the side of making 
a conservative overestimate of the stray light. The combined differential scattering and 
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diffraction can thus be expressed as 



di 



X 



A 




( 



A2 



) 



(7) 



where — a^max — a^min is the illuminated mirror width, 9^ is the non-specular 
angle measured in the plane of Figure |3] is a wing enhancement factor that takes 
account of deviations from a perfect diffraction aperture, a is the root-mean-squared 
(r.m.s.) micro-roughness of the mirror, and A is a horizontal coherence length for the 
mirror's non-ideal surface modulations. The first term in square brackets above takes 
account of diffraction due to the mirror's finite size, and the second term accounts for 
micro-roughness. 

Equation dTji assumes both that the micro-roughness is small compared to the wave- 
length {i.e., cr <S A) and that the non-specular angles are large compared to the angular 
size of the specular beam {i.e., 9^ ^ X/Dx). We denote the relative fraction of the 
diffracted solar-disk light on the mirror surface that makes it through the slit as cl///o. 
This fraction is given by integrating over all non-specular angles between 6'min and 
6'max, and 



Note that the integration over the other non-specular angle 9y is discussed in the Ap- 
pendix. 

5. Flux Normalization 

Next, it is necessary to compute the total flux that comes from a given location on the 
mirror and passes through a given part of the slit. We must take account of the full two- 
dimensional shape of the slit aperture; i.e., not just the width [w] but also the height [h] 
of the spatial element of interest (in the direction normal to the page in Figure[T]l. From 
the point of view of the mirror, the solid angle of the slit is given by ftsut = wh/ f"^. 
The flux of light from the extended corona is thus given by 



where /cor is the coronal intensity, in units of erg s~^ cm"^ sr~^ {i.e., integrated across 
the spectral line). We assume that /cor is constant over the slit aperture. 

The stray light, which began its life as the direct solar-disk intensity Iq, is attenu- 
ated by diffraction and non-specular scattering. Without this attenuation, the flux from 
the solar disk is given by /©il©, where is the solid angle of the Sun in the sky 
(essentially tt^q). However, because of the finite slit height h, only a fraction of this 
flux would make it through the slit. In the limiting case where the angle subtended by 
the slit-height direction (approximately h/ f) is small compared to the solar diameter 
{Sq = 29q), the piece of the solar disk that is "seen" through the slit is only a strip with 
solid angle ilstrip ~ It-Sq/ f. In the opposite limit of /i// > Sq, then f^strip = ^q- The 




(8) 




(9) 
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code used to compute the stray light uses an exact expression that bridges these hmiting 
cases, which for h/ f < Sq is given by 



(10) 



For all of the numerical calculations shown below, however, the approximation fistiip ~ 
HSq //is reasonably valid. Thus, the fully attenuated stray-light flux is 

The two fluxes defined above represent power per unit mirror surface area. The total 
power, in erg s^^, that passes through the slit for the two cases is given by integrating 
over the mirror surface: 

Poor = J Jdxdy F,,,{x) (12) 

/ dxdy Fstrayix) (13) 

where the integration over the y-coordinate (the mirror dimension out of the page in 
Figure [U is trivial because the fluxes are independent of y. The x-coordinate ranges 
from Xinin to Xmax- The y-coordinate ranges from y^nin = to j/max = 5 cm (the latter 
being the UVCS mirror height). 

Finally, the "calibrated" measurement for stray-light intensity can be constructed. 
Once the total power quantities defined above are known, the corresponding intensities 
can be computed simply by dividing either Pcor or Pstray by the product AJ7siit . The 
illuminated mirror area is given by 

^ — (ymax ymin)(*^max •^min) ■ (14) 

This process is unnecessary for the coronal intensity, since /cor was used in its initial 
definition [Equation (|9])]. However, the stray-light intensity [/stray] is computed in this 
way by dividing Pstray by ^Osiit directly. We do not consider any instrumental proper- 
ties beyond the entrance slit (e.g., efficiencies of the grating and detector) because they 
do not distinguish between coronal and stray light. 



6. Results 

We implemented the above constraints on ray geometry, external-occulter diffraction, 
non-specular scattering, and flux calibration into an IDL code that predicts the stray- 
light intensity [/stray] for various configurations of the UVCS instrument. A standard 
set of models for the H I Lya spectral line was constructed with the geometric param- 
eters given in the previous sections, as well as the following choices for the internal 
occulter position (xint = 1-5 mm), sUt width (w — 100 /im), and the angular resolution 
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of the spatial element of interest along the slit (10 arcseconds, or h — 36.4 /im on the 
slit). 

In the models presented here, the two main properties of the mirror that describe 
its diffraction and non-specular scattering are the r.m.s. micro-roughness [a] and the 
r.m.s. wavefront error [w]. As described in the Appendix, the wavefront error is a single 
parameter that characterizes the large-scale figure error of the mirror This parameter 
determines the angular width of the specular beam and the enhancement of the large- 
angle diffraction wing. The figure error can also be expressed as an effective mirror size 
-Dcff, which acts as a large-scale coherence length analogous to the micro-roughness 
coherence length [A]. Although no precise measurements of these quantities exist for 
the UVCS mirrors, there are several reasonable lower and upper bounds that can be 
applied. 

First, a "best-case" model can be constructed with the assumption of zero micro- 
roughness. For (7 = 0, the non-specular rays are solely the result of diffraction around 
the finite aperture of the primary mirror The lowest possible level of non-specular 
intensity is consistent with the diffraction-limited case of no figure error (w — 0). This 
limit is interesting because it provides a true lower bound on the stray light that cannot 
be improved upon without changing the instrument design. 

Second, we constructed a model with empirical upper limits on the micro-roughness 
and figure error. Scattering measurements were made with a helium-neon laser on the 
mirror eventually chosen for flight as the Lya channel primary ( |Leviton et al, 1998] ). 
These measurements were consistent with a value for the r.m.s. micro-roughness of 
approximately 20 A. We note that this determination is dependent on a particular nu- 
merical model for the scattering properties of mirrors, which has been described by 
ISaha et 011(119961) . Nonetheless, we use cr w 20 A as an approximate upper limit for 
the UVCS primary mirrors. Two independent determinations of the spatial coherence 
lengths of the mirror were found that agree with one another reasonably well: 

We used the lLeviton et a/. Id 1998] ) measurements of non-specular scattering to put a 
limit on the mirror's coherence length. The large-angle limit of Equation (|26] | can be 
used to show that 



Measurements of d///odr2 were taken at angles {9} corresponding to heliocentric 
distances between 0.95 and 11.3 Rq. The mean value of A for all of these measure- 
ments, derived from Equation ( fTSb , was found to be 1.36 cm. This results in a ratio 
[Dx/ K\ =5.3 for the full unvignetted mirror width = 7.2 cm. 

ii) Another determination of the mirror's performance can be found from the measured 
UVCS spatial resolution. The most definitive upper limits on spatial resolution were 
obtained from in-flight observations of the star p Leo, made in August 1996. For H I 
Lya, the measured full-width at half-maximum (FWHM) of the spatial resolution 
element was found to be 5.7 arcseconds. Although some fraction of this value may 
be attributable to other instrumental factors, it is useful to assign it all to the mirror's 
figure error to obtain an upper limit. Ideally, the FWHM should be converted into 
an angular half-width [9s] as defined in the Appendix. If the measured profile were 
a Gaussian beam, it would encompass 83.8% of the encircled intensity out to an 




SOLA: cranmer_stray0 9.tex; 21 January 2010; 16:37; p. 10 



UVCS Stray Light Contamination 




Figure 4. Comparison of modeled and observed stray-liglit suppression ratios /stray / Iq ■ Models are shown 
for the lower-limit case of no micro-roughness and diffraction-limited optics (lower solid curve) and an 
empirical upper limit on the mirror parameters (upper solid curve). The stray-light measurements (points) 
are described in the text. For comparison, the H I hya intensity measured by UVCS in a polar coronal hole 
is also shown (dashed curve). 



angular half-width of about 0.8 times the FWHM. Because this is so close to unity, 
and because its exact shape is unknown, we just use the FWHM itself to assume an 
upper-limit value of Og ~ 5.7 arcseconds. Thus, Equation ( |32] | gives Deff ~ 0.54 
cm. The stellar observations were taken at a heliocentric distance robs ~ 5.2 Rq at 
which the unvignetted mirror width was w 3.03 cm. The derived ratio Dx/Dos 
was found to be 5.6. 

Because the ratios I?x/A and D^/Dcs were found to be in good agreement with one 
another, we make the tentative assumption that one can treat A and Dctf roughly inter- 
changeably for this analysis. Thus, we adopt a mean ratio Dx/ A. = 5.45 and we use 
Figure|5]to give the corresponding wavefront error w/A w 0.23 and wing enhancement 
ratio w 1.78. 

Figure|4]shows the numerically computed /stray / /© ratios that correspond to the two 
limiting sets of parameters described above. The lower-limit set of mirror parameters 
(a = 0, Dx/A — 1) shows the effect of having only diffraction contribute to the 
stray light. For the empirical upper- limit parameters (cr = 20 A, Dx/A = 5.45) the 
micro-roughness contributed roughly half (45%) of the total stray light. The following 
parameterized fits were found to be good representations of the model curves in Figure 

il 

/ /stray \ ^ 3.14 X 10^« 

V ^0 Aower ([r//?0] - 1.121)4-56 ' ^ ^ 
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/stray \ ^ 1-30 X 10'" 

/upper ([r/i?o] - 1-076)4.73 ' ^ ^ 

and these fits apply mainly to heights between 1.45 and 5 i?©. 

The measured stray-light data points in Figure|4]were obtained from several sources. 
Pre-launch laboratory testing provided the upper limit value of 10^^ (filled circle) at 
r = 2.7 i?0 (jYsMetal, 1995| l. Observations of the faint Si III 1206 A emission Une - 
which has no coronal counterpart and is 100% stray light when measured by UVCS - 
were made in June 1996 (triangles) and scaled to the disk intensity of H I Lya to obtain 
the stray-light suppression ratio ( [Gardner ef a/., 1996) 1. Additional data points obtained 
from measurements of the C III 977 A line (squares) are also shown. Several of these 
measurements were discussed bv lKohl et a/. I (1 19991) . The remaining measurements are 
those at 1.64 Rq ( [Raymond et al.,2QQ3) and at 3.5 i?© ( [Suleiman etal, 1999) 1. In the 



latter two cases, only the total intensities were reported in the cited papers. These values 
were normaUzed to the solar disk and scaled to H I Lya using the relative line inten- 
sities of lVernazza and Reevesl ( [T978) l and the time-dependent irradiance measurements 
of [Woods et al.\ (I2000 i The error bars on the C III points were estimated from the 
assumption of an overall 50% uncertainty, which can be split up evenly between 25% 
uncertainties in the numerator [/stray] and the denominator [/©]. 

Figure H) does not show data for distances larger than robs ~ ARq. At these large 
observation heights, the measured stray light is dominated by a roughly constant back- 
ground. This background is likely to be the result of several effects, including the 
dust-scattered F-corona, interstellar line emission (for H I Lya and other lines of neutral 
species), and possibly other sources of instrumental stray light that are negligible at 
lower heights. In visible-light Lyot coronagraphs, a constant stray-light background of 
this kind is often the result of dust particles on the primary mirror ( [Nelson et al, 2008) ). 
However, for the externally-occulted design of UVCS - in which the relatively small 
primary mirror is never exposed to direct sunlight - the effects of dust have been 
considered negligible at both visible and ultraviolet wavelengths ( [Romoli ef a/., 1993) 
[Kohl et al, 1995 1). The vast majority of UVCS coronal measurements have been made 
within the range of heights shown in Figure |4) 

For any given observation height, the intensity of the extended corona can vary by 
several orders of magnitude. Figure|4)shows the mean radial dependence of the intensity 
of H I Lya in polar coronal holes measured during the 1996 - 1997 minimum of solar 
activity ( [Cranmer ef a/., 1999) 1. Coronal holes are the lowest-intensity structures in the 
extended corona, and thus they represent the most stringent requirements on stray-light 
suppression. The coronal intensities were normalized by dividing by a disk intensity 
typical for the time period of these solar-minimum observations, /© = 5.3 x 10^^ 
photons s^^ cm^^ sr^^ {e.g., I Woods et al. \ 120001) . Even the emission from coronal 
holes is significantly larger than the modeled and measured stray-light levels for H I 
Lya. Other spectral lines tend to exhibit larger ratios of /cor//© than H I Lya, so the 
stray light is less of a contaminant. 

In order to scale the above results to other spectral lines, the overall wavelength 
dependence of /stray must be calculated. The diffraction from the external occulter 
depends on A because of the definition of v in Equation ([3]). Also, the total rate of 
diffraction and scattering through the slit [Equation ()7)l] scales differently with A de- 
pending on whether the non-specular stray light is dominated by finite mirror-size 
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diffraction or by micro-roughness. At a standard observation height of 2 Rq, we com- 
puted /stray for a range of wavelengths between half and twice that of H I Lya. The 
lower-limit case of u = cj = gives rise to an approximate wavelength depen- 
dence /stray oc A^'^^. The empirical upper-limit choice of mirror parameters gives 
/stray oc A° '^'*. These exponents increase slowly as a function of increasing height. 
For the bright O VI 1032, 1037 A doublet, there is only a 15% relative separation in 
wavelength from 1216 A. Since the uncertainties in the model parameters exceed this 
level at any given height, it is reasonable to just use the modeled value of /stray that was 
computed for H I Lya. 



7. Conclusions 

The primary aim of this paper has been to construct a model of the ultraviolet stray- 
light properties of the UVCS instrument on SOHO. The model has essentially no freely 
adjustable parameters, with the possible exception of the parameters a and oj that define 
the mirror imperfections. Estimates of these parameters were obtained from pre-launch 
laboratory tests and in-flight measurements. A completely separate set of existing mea- 
surements of the stray-light suppression ratio for H I Lya was used to test the model. As 
can be seen in Figure|4l the modeled stray light agrees well with these measurements. 

It is interesting to investigate how much of an improvement in stray-light suppres- 
sion can be made by various changes to the externally occulted coronagraph design 
exemplified by UVCS. Equation Q indicates that even if the mirror were made to be 
perfectly smooth (a ~ 0) and ideally diffraction-limited (W — 1, = A), there 
would still be a finite level of non-specular stray light due to diffraction by the finite- 
sized mirror. Reductions in the amount of non-specular radiation can then be achieved 
by increasing the mirror size [D^] or by increasing the magnitudes of the non-specular 
angles [9x] that pass through the slit. 

Designs for next-generation coronagraph spectrometers have also achieved improve- 
ments in performance by increasing the distance [//oxo] between the mirror and external 
occulter ( |Kohl et al, 2006[|Kohl et al, 2008| l. An increase in Hcxo results in a smaller 
value of //diff [Equation (|2|i], which leads to smaller non-specular angles [9x]- Taken 
by itself, this would give rise to an increase in the non-specular intensity that scatters 
and diffracts from the mirror into the slit. However, this effect is more than offset by 
the larger unvignetted mirror area that becomes available when //cxo is increased (i.e., 
an increase in X2). The larger mirror area gives rise to a "deeper shadow" under the 
external occulter, which is equivalent to a larger extent of v in Figure |2] This results in 
a lower overall amount of diffracted light striking the mirror. Also, the larger available 
mirror area allows more coronal rays to be collected, thus increasing the overall instru- 
mental sensitivity. These advantages, together with improvements in mirror reflectivity 
and detector efficiency, can give rise to several orders of magnitude improvement in 
sensitivity compared to UVCS/SOHO. 
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Appendix 

A. Derivation of Non-Specular Intensity Profiles 

Here we describe the angular dependence of non-specular diffraction and scattering that 
arise from reflection (at nearly normal incidence) from a mirror of finite size and with 
a non-ideal surface. The fraction of the incident intensity /q scattered into a differential 
solid angle [dl7] is defined here as 

where the scattering wavevector q is the vector difference between the incoming and 
reflected photon wavevectors k. Thus, if the mirror surface is assumed to lie in the x, y 
plane, then {qx,qy,qz) ~ "^T^i^x, Sy, 2)/A. The goal of this Appendix is to show how S 
depends explicitly on the non-specular angles 9x and 9y. The specularly reflected ray is 
defined as 9x — 0y — 0- 

In the optics literature, the differential scattering fraction defined above is often 
described using other terms. A frequently used variant is the bi-directional reflectance 
distribution function (BRDF), which is defined as 

BRDF = , (19) 

cos H lodil 

and where 8 is the scattering angle measured from the local normal to the mirror surface 
( [Bennett and Mattsson, 1 999) 1 . Another commonly used term is the power spectral den- 
sity (PSD), which is a measure of the angular distribution of fluctuations on the mirror 
surface itself. The PSD can be expressed in similar units as the angular distribution of 
reflected power, and 

\4 J J 

The dimensionless quantity TZ is the overall reflectivity of the mirror, and it is often 
called an "optical factor" that depends on wavelength, the angles of incidence and re- 
flection, and the degree of polarization of the incident and reflected beams ( [Bennett, 2003[ ).l 
In order to separate out the problem of mirror reflectivity from the problems of non- 
specular scattering and diffraction, the analysis below makes the assumption that TZ = 
1. 

The angular intensity profile can be derived by applying the Born approximation 
from scattering theory (e.g. JSinha et aZ.i ri988 ; Church and Takacs 1993 ; Gu llikson e t al. | 
I1997I I. If the reflectivity of the mirror is assumed to be perfect, the Born approximation 
can be expressed as 

^(q) = ^JdxdyJ dx' dy' e-'^-^''"''^ . (21) 

Both sets of spatial integrals are taken over the mirror area A in the x, y plane. The 
position vector r depends not only on the horizontal position, but also on the spatially 
varying mirror height z{x, y). 
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In the UVCS instrument, the most important orientation for the non-specular diffrac- 
tion and scattering is the Ox direction {i.e., the left-right horizontal direction in Figures 
[U and O. Because light from a given point in the y-direction (on the sky or on the 
external occulter) fills the mirror in the y-direction, there ends up being a wide range 
of 6y scattering angles that go through the slit. Thus, we first solve Equation (1211 1 and 
then integrate over all 6y angles to obtain the differential scattering rate d//(/od6'2.). 
Sections A.l and A.2 give independent calculations of this quantity in two limiting 
regimes of the spatial scales of mirror imperfections. Equation ^ above estimates their 
combined effect as the sum of the two rates. 



A. 1 . Non-Specular Scattering due to Micro-Roughness 



When only small-scale statistical fluctuations in the mirror surface are taken into ac- 
count, it is possible to describe the mirror height with spatially averaged moments. The 
mean height is defined as 



(z) = ^ y dxdyz{x,y) = , 



(22) 



and the rm.s. surface roughness is given by (z^)^/^ = a. As defined here, a contains 
no information about the horizontal distribution of microscopic structures on the mirror 
surface. Typically, laboratory measurements of <j are dependent on the range of horizon- 
tal scales that were probed. Many real surfaces have a "fractal" distribution of surface 
fluctuations that are ill-described by a single a. To begin the process of characterizing 
the horizontal scales, though, we define an autocorrelation function 

C{s,t) — ~^ J <ix dy z{x,y) z{x — s,y — t) (23) 

and we assume the surface is isotropic on small enough scales such that C is a function 
of only r = (s^ +t2)i/2. 

With the above definitions, the scattering integrals in Equation (1211 1 can be shown to 
reduce to 

SiO) ^^Jdrre.^ [-'-^ ^ - } Jo (^) (24) 

where Jq is the zeroth-order Bessel function of the first kind, and 6 — {Ox+Oy)^^^ is the 
full magnitude of the scattering angle away from the specular beam. In this section we 
assume the integration over r can be taken from zero to infinity, which ignores the finite 
size of the mirror This is equivalent to assuming that the patterns of micro-roughness 
are small in comparison to the mirror size. We can make use of a commonly used (and 
empirically derived) parameterization for the autocorrelation function, 

C{r) ^ a^exp[~{2TTr/Af''] , (25) 

in which A is a horizontal correlation wavelength (or coherence length scale) for the 
mirror and /i is a texture parameter that describes the shapes of the surface features. 
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Small values {h <^ 1) indicate extremely jagged surfaces and large values {h « 1) de- 
scribe smooth crests and troughs. Sometimes the autocorrelation function is expressed 
as a sum of several terms each similar in form to Equation dZSl ). but with different 
values of h ( |Elson et al, 1993) 1. In that case, the rm.s. surface roughness a in each 
term can refer to different ranges of spatial scale. The special value h = 1/2 is often 
identified as a "conventional surface finish," and we will restrict ourselves to this value 
for simplicity (see ISinha etaU [T9881 IChurch and Takacsl [T9931 IFineschi et aU IT994] |. 
Using this value, combined with the approximation of smaU-scale surface roughness 
(cr <C A), the intensity scattering profile becomes 

~ 2^ + A4[l + (0A/A)2]3/2 ■ (26) 

The first term, proportional to a Dirac delta function, describes the ideal specular beam. 
The second term is the scattering due to micro-roughness, and its integral over all solid 
angles gives the so-called total integrated scatter (TIS), or 



TIS = 27r / deeS(e) = ^^^^ . (27) 



167rV2 

Noting that 9^ = + assume 9x and integrate over 9y to obtain 

dl levrcr^A 



(28) 



(29) 



Iod9x A3[1 + (0,A/A)2] ■ 
For large angles, far from the specular beam (i.e., 9^ ^ A/A) this reduces to 

dl _ 167ro-2 

which is used in Equation (|7]i above. 

A. 2. Non-Specular Diffraction due to Finite Mirror Size and Aberrations 

For an illuminated mirror that reflects a finite-sized (rectangular) incident beam of 
light, diffraction will broaden the reflected beam and give rise to power at non-specular 
angles. This situation is similar to the diffraction of parallel rays that pass through a 
rectangular aperture (see, e.g.. lSchroeder 2000, H echtl l20021 l. 

The simplest case of a mirror with an ideal surface can be treated by assuming z — Q 
in Equation dTTT i and taking the integration limits to be finite in both x and y. This yields 

DxDy I sin \ / sinw. 



where = (x^ax - a^min) and Dy = (j/,„ax - ymin) are the dimensions of the 
illuminated part of the mirror The dimensionless angle coordinates are defined as Ui = 
TTDi9i/\, for i — x,y. Equation ( l30l l is the well-known diffraction pattern that drops 
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off as u ^ along the normal axes of the rectangular aperture, and as w along the 
diagonals. As above, we integrate over 9y to obtain 



dl Dt ( sinu 



/ode, \ \ ) 2^2^,02 



(31) 



The latter approximation above is taken in the limiting case of large-angle scattering 
(itx 3> 1), where the mean behavior of the function (sin w/m)2 can be approximated as 
l/(2«2). 

For diffraction off a mirror containing large-scale aberrations, we must again apply 
a variable mirror height z(x, y) in Equation (|2TI) . This is traditionally done by inserting 
a phase factor that takes account of the optical-path differences between the aber- 
rated wavefronts and the ideal, non-aberrated wavefronts. IBorn and Wolfl (1 19991 1 and 
ISchroederl (j2000 ) defined these phase factors in terms of Zernike polynomials for var- 
ious types of common mirror deformation patterns. Aberrations change the diffraction 
profile defined above in two distinct ways: 

They enhance the non-specular scattering in the far wings of the profile (see, e.g.. 
Figures 10.12 to 10.17 of ' Schroederl I2000 i Here, we take this into account by 
multiplying Equation OTT i by a wing enhancement factor W . 

ii) They broaden the specular beam, which for a finite mirror or aperture is defined ac- 
cording to how far from the ideal specular ray one must go to encompass a specified 
level of "encircled energy." The broadening of the specular beam is related closely 
to how far the actual mirror is from the diffraction limit. 

In order to quantify these effects, we calculated a range of diffraction profiles for various 
types of aberration. Kohl et aZ. I (1 1995 *) noted that spherical aberrations tended to be the 
dominant type for the UVCS mirrors. The remainder of this section concentrates only 
on that type. The standard measure of the magnitude of an aberration is the r.m.s. wave- 
front error uj, which is usually expressed as a dimensionless number of wavelengths 
w/A. The relative departures from the perfect (non-aberrated) case were computed for 
a circular mirror, but they should be comparable for mirrors of other shapes. 

Figure |5] shows both the broadening of the specular beam and the enhancement of 
the power-law scattering wing as a function of the wavefront error. The width of the 
specular beam was computed by determining the angular radius [9s] at which the encir- 
cled energy was 0.838 times the total energy distributed over all angles. This fraction is 
equivalent to the light enclosed by the the classical 'Airy disk" (see H echt 2002). For 
a diffraction-hmited circular mirror without aberrations, 9s — 1.22X/D. A mirror with 
aberrations has a larger value of 9s and thus corresponds to a smaller effective diameter 
Dcs, which would have given rise a specular beam of identical width in the absence 
of aberrations. Once Og has been computed for a given magnitude of the spherical 
aberration, it is then possible to solve for the ratio 

D 9, 

> 1 . (32) 



Doff 1.22 A/Z) 

The quantity Doff characterizes the horizontal spatial scales of large wavefront de- 
viations on the mirror surface, just as A characterizes the coherence length of the 
micro-scale wavefront deviations. Interestingly, laboratory and in-flight measurements 
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Figure 5. Effects of spherical aberrations on the ratio of the mirror diameter D to the effective diameter D^ff 
(solid curve) and on the relative wing enliancement W (dashed curve). 

for the UVCS Lya primary mirror yielded similar values for the ratios Dx/Dcs and 
Dx/A (see Section 6). 

We computed the magnitude of the enhanced diffraction wing by fitting the angular 
dependence of S{9) with a power-law (outside the specular beam). The ratio W was 
found by taking the ratio of the numerically computed wing intensity to that produced 
by a mirror with no aberrations. The dependence of this ratio is also shown in Figure 
|5] Thus, if any one of the three quantities (oj/X, D/Dcs, or W) is specified for a given 
mirror, Figure|5]can be used to estimate the other two. 
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